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ABSTRACT 

Solar eruptions, particularly coronal mass ejections (CMEs) and extreme-ultraviolet (EUV) jets, have 
rarely been investigated with spectroscopic observations. We analyze several data sets obtained by the 
EUV Imaging Spectrometer onboard Hinode and find various types of flows during CMEs and jet erup- 
tions. CME-induced dimming regions are found to be characterized by significant blueshift and enhanced 
line width by using a single Gaussian fit. While a red-blue (RB) asymmetry analysis and a RB -guided 
double Gaussian fit of the coronal line profiles indicate that these are likely caused by the superposition of 
a strong background emission component and a relatively weak (^10%) high-speed (~100 km s _1 ) up- 
flow component. This finding suggests that the outflow velocity in the dimming region is probably of the 
order of 100 km s _1 , not ~20 km s _1 as reported previously. Density and temperature diagnostics of the 
dimming region suggest that dimming is primarily an effect of density decrease rather than temperature 
change. The mass losses in dimming regions as estimated from different methods are roughly consistent 
with each other and they are 20%-60% of the masses of the associated CMEs. With the guide of RB 
asymmetry analysis, we also find several temperature-dependent outflows (speed increases with tempera- 
ture) immediately outside the (deepest) dimming region. These outflows may be evaporation flows which 
are caused by the enhanced thermal conduction or nonthermal electron beams along reconnecting field 
lines, or induced by the interaction between the opened field lines in the dimming region and the closed 
loops in the surrounding plage region. In an erupted CME loop and an EUV jet, profiles of emission 
lines formed at coronal and transition region temperatures are found to exhibit two well- separated compo- 
nents, an almost stationary component accounting for the background emission and a highly blueshifted 
(~200 km s _1 ) component representing emission from the erupting material. The two components can 
easily be decomposed through a double Gaussian fit and we can diagnose the electron density, tempera- 
ture and mass of the ejecta. Combining the speed of the blueshifted component and the projected speed 
of the erupting material derived from simultaneous imaging observations, we can calculate the real speed 
of the ejecta. 

Subject headings: Sun: coronal mass ejections (CMEs) — Sun: flares — Sun: corona — line: profiles — solar wind 



1. Introduction 

Coronal mass ejections (CMEs) are large-scale 
solar eruptions and earth-directed CME s are often 
sourc es of strong g eomagnetic s t orms (e.g., Gosling et al.1 
19911: IWang et alJ I2002L 120061: Izhang & Lowl 120051: 



Feng et al 1 120091: iLiu etaD 1201 Ob. R ece nt statistical 



studie s of lReinard & Biesecker^ 2008) and lBewsher et al 



(2008) have shown that more than 50% of the frontside 
CMEs are associated with coronal dimmings (or 
transient coronal holes), which are characterized 
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(EUV) and soft X-rays (e.g., 


Rust &Hildnerll 19761: 
Thompson etal.ll998: 


GoDalswamv & Hanaokall998 


Zarro et al.1 Il999t Izhou et al. 


2003: beTomaetal. 


20051: IMcIntosh et alJ 120071 Miklenic et al.1 2011). 



Dimmings may mark l ocations of the fo otpoints of 
ejected flux ropes (e.g.. IWebb et al.l 12000) or formed 
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by reconnection between the erupting fie ld and the 



surrounding magnetic stru ctures (e.g., lAttrill et al 
2007; Mandri ni etaDl2007l) . There are basically two 
types of dimmings: small-scale dimmings associated 
with t he two ends of a pre-CME sigmoid structure 
(e.g.. ISterling& Hudson! 1 19971: [ Hudson et alJ I l998[ 



Zarro et alJI 19991: [Webb et alJbOOd: Ijiang et al J 120031: 
Cheng et alJEoToh and global-scale dimmings which 



are often i mmediately proceeded by global "EUV 
waves" (e.g jThompson et al.l20o5 lAttrill et al.l2007l) . 

Jet-like phenomena are small-scale solar erup- 
tions and they are often observed in X-ray, EUV, and 
white light. Most jets a re associated with small flares 
dMadjarska et al.lEo07|) . EUV jets are characterized 
by nearly collimated high-speed motions of plasma at 
coronal and transition region (TR) temperatures (e.g. , 
Alexander ^et^l 19991; K^T^bOOdlUu et al 11201 ll: 
Shen et al JI20TH l2012t Isirivastava & MurawskilEoTTl) . 
Studies have shown that EUV jets a nd X-ray jets are 
closely associated with each other (IKim et al.l 120071 



Chifor et al. 2008b; He et alJ2 010b: Y ang et al.ll201ll) . 



Recent observa tions by the Atmos pheric Imaging As- 
sembly (AIA, Lemenet "all l201lh onboard the Solar 
Dynamics Observatory (SDO) have revealed that fine- 
scale EUV jets (high-speed outflows) are ubiquitous 
on the Sun dDe Pontieu et al.ll201ll : iTian et alJEoilbl 



Yang etal.ll2011h . 



Kinematics associated with CMEs and EUV jets are 
usually studied through coronagraph and broadband 
observations. The high cadence and large field of view 
(FOV) of these imaging observations have greatly en- 
hanced our understanding of these solar eruptions. 
However, imaging observations only allow us to study 
the plane of sky (POS) component of the kinematics, 
which is usually a good approximation of the full kine- 
matics only for limb events. For earth-directed erup- 
tions, especially halo-CMEs which are the cause of 
most strong geomagnetic storms, imaging instruments 
placed close to the Sun-Earth line often fail to observe 
their initial or complete evolution. Spectroscopic ob- 
servations, on the other hand, can provide information 
on the plasma motions in the line of sight (LOS) di- 
rection and thus are critical for us to understand the 
kinematics of earth-directed eruptions. For both disk 
and limb eruptions, their three-dimensional (3-D) evo- 
lution can in principle be revealed through simultane- 
ous imaging and spectroscopic observations. In addi- 
tion, spectra of different emission lines can be used 
to diagnose plasma properties such as electron density 
and temperature. Spectroscopic data can also be used 



to estimate the mass of the erupted material and mass 
loss in the dimming region. 

So far there are only a few spectroscopic investi- 
gations of CMEs, dimmings, and EUV jets in the lit- 
erature. Using observati ons by the Corona l Diagnos- 
tic Spectrometer (CDS. lHarrison et al.ll 19951) onboard 
the Solar and Heliosp heric Observatory (SOHO), 
Harra & Sterling (2001) reported significant blue shift 
of emission lines formed at coronal and TR tem- 
peratures in dimming regions. This result has been 
confirmed by recent high-resolution o bservations of 
the E UV Imaging Spect rometer (EIS. ICulhane et al 
20071) onboard Hinode faarra et aD 120071 Ijin et al 



2009[lAttrill et al.l2010l:IChen et al.l201 Pal IHarra et al 
201 lah . Blue shift was also f ound in footpoint regions 
of small-scale erupted loops dHe etal.ll2010bl) . How- 
ever, a preliminary study of iMcIntosh et al.1 ( 2010|) 
suggests that some line profiles in the dimming re- 
gions are asymmetric, with a weak enhancement in the 
blue wings. EIS observations have also revealed an 
obvious increase of the line broadening in dimming 
regions, which was interpreted as a growth of Alfven 
wave amplitude or inhomogeneiti e s of flow velocities 
along the LOS (iMcIntoshl 120091: IChen et al.l EoiOal 
Dolla & Zhukovl l201ll) . The presence of asymmet- 
ric line profiles suggests that there are probably two 
emission components and that a single Gaussian fit 
may not reveal the real physics in dim ming regions 
dMcIntosh et al .11201 Pi bolla & Zhukovll201ll) . 

The outflow speed derived from a single Gaussian 
fit is roughly in the range of 10-40 km s _1 and usually 
it does not change significantly for coron al emission 



lines formed at different temperatures (IHarra et al 



20071: Ijin et all 120091: lAttrill et al.l l2010l IChen et al 



2010a|) . However, llmada et aD ( 2007|) reported a 
temperature-dependent outflow in the dimming region 
following a CME. The speed of the flow increases 
from -10 km s" 1 at log (77K)=4.9 to -150 km s" 1 at 
log (77K)=6.3. One-dimensional modeling effort has 
been tak en to reconstruct this temperature-dependent 
outflow (llmada et al.ll201lk . 

Line splitting is usually associated with a very high- 
speed (—200 km s _1 or large r) plasma motion Usin g 
CDS and EIS ob servat ions, IHarra & Sterlingl d2003l) . 
lAsai et al.l d2008h and iLi & Dingl d2012l) found sig- 
natures of line splitting indicative of plasma ejection 
at a speed of —250 km s _1 during CMEs or fila- 
ment eruptions. Spectra obtained by the Solar Ul- 
traviolet Measurem ents of Emitted Radiation Spec- 
trograph fSUMER. Iwilhelm et al.ll995l:lLemaire et al 
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1997k onboard SOHO have revealed signatures of line 
splitting associated with the expanding X -ray plasma 
in a flare/CME event dlnnes et al.ll200l|) . Line split- 
ting or obviously blueshifted components have also 
been fou nd in spectra of EUV jets in coronal holes 
and ARs (w ilhelm et al 1 120021 iMadiarska et"aDl2QQ7l: 
Kamio et al .1120071 12009l1chif or et al.ll2008ah . 

There have been a few investigati ons of the plasma 
proper ties of dimmings and EUV jets. lHarrison & Lyons 
(l2000h and lHarrison et al.ld2003h used the Si x 347.40A & 356. 
line pairs to diagnose the electron density and found 
that it decreased as dimming occurred. Using some as- 
sumptions of the emitting volume and the distribution 
of the amount of material at different temperatures, 
they also made an effort to estimate the mass loss in 
the dimming region and found that it is of the same 
order as the mass of the associated CMEs. Taking 
values of the formation heights of different emission 
lines an d the densities from static solar atmosphere 
models, Ijin et al.1 (2009) also developed a method to 
estimate the mass losses in dimming regions associ- 
ated with two events during 2006 Dec 1 3-15. Using 



emission often consists of at least two components so 
that previous results based on a single Gaussian fit may 
need to be reconsidered. We apply the recently im- 
proved techniques of red-blue (RB) asym metry anal 



ysis a nd RB -guided double Gaussian fit (ITian et al 



the Fe x n 186.88A & 195. 12A line pair. lChifor etal 



(2008a) measured electron densities higher than log 
(7V e /cm -3 )=ll for an EUV jet. However, they only 
simply summed up the spectral intensities in the wave- 
length windows of the lines and could not separate the 
blueshifted component from the background emission 
component. 

Besides dimmings and ejecta associated with CMEs 
and EUV jets, other solar eruption related phe- 
nomena such as flare induced chromospheric eva p- 
oration dTeriaca et al.1 120031: iMilligan et al.1 12006a 



Millisan & Dennis 2009; Millisan 2011; Chen et al 
2010b; Watanabe et al.l201(AlLi&Ding 201 1; Gra ham et al 



20111: iLiu et al. 201 JJ), flare-related m agnetic recon- 
nection dWang et alJl2007l: lHara et al.ll201 lk, filamen t 
oscillations dChen et al.l f2008: Bocchialini et al.ll201ll) 



and coronal waves (jHarra et al.l 1201 lbc IChen et al.l 
201 lc IVeronig et al .11201 1|) have also been investigated 
through EUV spectroscopic observations. As the ap- 
proaching of the new solar maximum, there is no doubt 
that more spectroscopic observations will be employed 
to study solar eruptions since the high-resolution EIS 
instrument is still in good condition and the Interface 
Region Imaging Spectrograph (IRIS) is expected to be 
launched in 2012. 

In this paper we analyze several data sets obtained 
by EIS during CME eruptions and EUV jets. The 
shapes of the EIS spectral line profiles suggest that the 



l2011ch . which we used previously to study properties 
of the high-speed outflows in non-eruptive active re- 
gions (ARs), to the spectra acquired during solar erup- 
tions. We find various types of flows and discuss pos- 
jsjR[e mechanisms to produce these flows. We also di- 
agnose the density, temperature and mass loss (mass) 
of the dimming region as well as the ejected material. 
Our analyses demonstrate that EUV spectroscopic ob- 
servations can provide a lot of valuable information on 
solar eruptions. 

2. Observations, Single Gaussian fit, and RB 
asymmetry analysis 

Table [T] lists some of the observation details of the 
six events we analyzed. The class and peak time of the 
associated flare are also listed for each event. There 
were many fast repetitive rasters (with a scanning ca- 
dence of ~6 minutes) for Events 4&5 and we only 
present results for several of them in this paper. 

The SSW routine eis .prep. pro was applied to cor- 
rect and calibrate the EIS data. This includes CCD 
pedestal and dark current subtraction, cosmic ray re- 
moval, warm and hot pixels identification, absolute 
calibration, error estimation, and so on. The effects 
of slit tilt and orbital variation (thermal drift) were es- 
timated by using the SSW routine eis. wave .corr.pro 
and removed from the data. After that, a running av- 
erage over 3 pixels along the slit was applied to the 
spectra to improv e the signal to noise ratio. Note that 
Tian et al.1 ( 2011c|) used the median values of the mea- 
surement errors when averaging profiles over several 
pixels. In this paper we regard these line profiles as in- 
dependent measurements of the profile at a single pixel 
and use the uncertainty propagation theory to calculate 
the measurement errors for the averaged profile. This 
usually leads to smaller values of the errors, which is 
reasonable since the spatial average improves the sig- 
nal to noise ratio. 

As a common practice, a single Gaussian fit was 
applied to each spectrum. The line peak intensity, 
Doppler shift and line width can thus be derived. We 
assume zero shift of the profile averaged over each ob- 
servation region. We have to mention that the line 
width can be expressed in different formats and dif- 
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Peak Intensity (arbitrary unit) 



Velocity (km/s) 



Width (km/s) 



RB Asymmetry [70-130 km/s] 



2.48 3.01 3.55 4.09 4.63 5.16 5.70 - 42 -28 -14 14 28 4 2 44 50 56 62 68 74 80 -0. 21 -0.14 -0.07 0.00 0.07 0.14 0.2 1 
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Fig. 1. — Spatial distributions of the peak intensity, Doppler velocity and exponential width derived from the single 
Gaussian fit, and the average RB p asymmetry in the velocity interval of 70-130 km s _1 for Fe XIII 202. 04 A in the 2006 
Dec 14-15 observations. The beginning and ending time (hour:minute) of each scan is indicated in the intensity image. 
The pre-emption conditions are shown in the first row. The square in each panel of the second row marks locations 
where profiles are averaged and presented in Figure [71 The red contours shown in the map of RB p asymmetry for the 
01:15-03:30 scan outline locations where the RB p asymmetry (70-130 km s _1 ) is smaller than -0.03 and the signal to 
noise ratio of the profile is larger than 8. 
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Table 1 : EIS Observations Used in This Study 



Obs. 
ID 


Scanning Period 


Exposure 

timp ( 


Slit 


Flare Class & Peak Time 


Comment 


1 


2006 Dec 14 15:11-16:01 
2006 Dec 14 19:20-21:35 
2006 Dec 15 01:15-03:30 
2006 Dec 15 04:10-06:25 
2006 Dec 15 10:29-11:19 


10 
30 
30 
30 
10 


1" 


X1.5,Dec 14 22:15 


CME & Dimming 


2 


2007 May 19 09:42-10:30 
2007 May 19 11:41-15:23 


10 
40 


1" 


B9.5,May 19 13:02 


CME & Dimming 


3 


2006 Dec 12 19:07-23:46 
2006 Dec 13 01:12-05:41 


30 
30 


1" 


X3.4, Dec 13 02:40 


CME & Dimming 


4 


2011 Jun 21 02:11-05:18 


9 


2" 


C7.7, Jun 21 03:25 


CME & Dimming 


5 


2011 Feb 14 19:13-20:06 


8 


2" 


C6.6, Feb 14 19:30 


CME 


6 


2007 Jun 5 01:51-02:04 
2007 Jun 5 04:16-04:29 


5 
5 


2" 


CI. 2, Jun 5 04:23 


EUVjet 



Peak Intensity (arbitrary unit) Velocity (km/s) Width (km/s) RB Asymmetry [70-130 km/s] 

2.48 3.01 3.55 4.09 4.63 5.16 5.70 -42 -28 -14 14 28 42 44 50 56 62 68 74 80 -0.21 -0.14 -0.07 0.00 0.07 0.14 0.21 




150 



X (arcsec) X (arcsec) X (arcsec) X (arcsec) 

Fig. 2. — Same as Figure [T]but for the 2007 May 19 observations. The rectangular regions 1 & 2 mark the locations 
where the Fe XIII 202. 04 A profiles are averaged and presented in Figure 



5 



Peak Intensity (arbitrary unit) 



Velocity (km/s) 



Width (km/s) 



RB Asymmetry [70-130 km/s] 





-50 


sec) 


-100 


(arc 


-150 


> 






-200 




-250 



D7-23:46 





-50 


o 

<D 
CO 


-100 


o 






-150 


> 






-200 




-250 




100 200 300 400 



01:12-05:41 








| 









100 200 300 400 
X (arcsec) 



100 200 300 400 
X (arcsec) 



100 200 300 400 
X (arcsec) 



100 200 300 400 
X (arcsec) 



Fig. 3. — Same as Figure \T\ but for the 2006 Dec 12-13 observations. The rectangular regions 1-4 mark the locations 
where line profiles are averaged and presented in Figures [T0lSfT4l 



ferent names are assigned t o different formats (e.g., 
Chae et al.lfl998l: |Peterll2Q10h . A Gaussian line profile 
can be expressed as: 



I(v) = I p exp(- 



1 (v- 
' 2 



(1) 



where v, I p and are the wavelength vector (con- 
verted into velocity thro ugh Dopple r effec t), peak in- 
tensity and line center. IChae et al.1 (11998b defined a 
as Gaussian width. While the Gaussian width men- 



tioned by |Pe^_^^js V2cr. In lTian et al.1 d2011cl) 
we followed [Peted d2010l) and used both the names of 
Gaussian width and 1/e width for To avoid con- 
fusion, in the following we use t he name e xponential 
width (or 1/e width, also used by lPeteri (|2010|) ) instead 
of Gaussian width for \/2a. 

The techni que of RB asymmetry an alysis was first 
introduced by iDe Pontieu et al.1 (|2009) and it is based 
on a comparison of the two wings of the line profile 
at same velocity ranges. The line profile was first in- 
terpolated to a spectral resolution ten times greater 
than the original one, then the blue wing emission in- 
tegrated over a narrow spectral range was subtracted 
from that over the same range in the red wing. The 
range of integration was then sequentially stepped 
outward from the line centroid to build an RB asym- 
metry profile (simply RB profile). In our previous 
work (IDe Pontieu et al. 2009: iDe Pontieu & Mcintosh! 
2010; De Pontieu et al. 2011; Mcintosh & De Pontieu 



d2011ch named this method RB $ and they further de- 
veloped two other methods RB p and RB p, , which are 
basically the same as RB# except the determination 
of the line centroid. For RBp, the spectral position 
corresponding to the peak intensity is used as the line 
centroid and the resulting RB profile is normalized to 
the peak intensity. For RBp,, the line center of the 
primary component, which is derived from the RBp- 
guided double Gaussian fit, is used as the line centroid 
and the resulting RB profile is normalized to the peak 
int ensity of the prim ary component. As pointed out 
bv lTianetal.1 d2011cl) . the RB p technique can resolve 
the blueshifted secondary component more accurately 
as compared to the originally defined RB# technique. 
Thus, here we apply the newly developed RBp and 
RBp, techniques, as well a s the RBp -guided double 
Gaussian fit (for details see iTian et al .11201 lcj) . to the 
data in this paper. 

Figures [T]|6] show the spatial distributions of the 
peak intensity, velocity and exponential width derived 
from the single Gaussian fit, and the average RBp 
asymmetry in the velocity interval of 70-130 km s _1 
for Fe XIII 202.04A or Fe XII 195. 12A in the observa- 
tions of six events. For the RBp asymmetry a neg- 
ative/positive value indicates an enhancement of the 
blue/red wing. We prefer to use the Fe XIII 202.04A 
line to detect asymmetry since there is no identified 
blends in this strong line, although the pervasive pres- 
ence of very weak redward asymmetries outside the 



2009aE;lMc_Intosh et al.|2oTT|;|Tian et al.|2011a|;|Martmez-SykqfMfeflyg regions (in Figures [IE]) might suggest an 



201 lh ~ we used the single Gaussian fit to determine the 
line centroid and applied this technique to spectra in 
coronal holes, quiet Sun, and quiet ARs. ITian et af 



unidentified weak blend at the red wing of the line pro- 
file. However, in Figures |4][6] we present results for the 
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Fig. 4. — Evolution of AIA 193 A intensity and EIS Fe XII 195. 12A line parameters (peak intensity, velocity and width 
derived from the single Gaussian fit, and the average RB p asymmetry in the velocity interval of 70-130 km s _1 ) in the 
201 1 Jun 21 observations. The time of the AIA observation and the beginning time of each EIS scan are indicated in 
the corresponding intensity images. The rectangular region marks the locations where line profiles are averaged and 
presented in Figure[TT] The size of the FOV is about 175" x 152". A movie (m4.mpeg) showing the evolution of AIA 
171 A is available online. 



■Vii ill l iiiiPih i'i-hi HnuniiF, V 1 1 ; .1 h r |itti! 



■Vl-i(irf4| l-N Jui h nwr u i h [ 7 fl i in hm.ii 




Fig. 5. — Same as Figure [4] but for the 2011 Feb 14 observations. The asterisk marks the center of five pixels along 
the slit where line profiles are averaged and presented in Figure [121 The size of the FOV is about 175" x 160". For 
illustration some bad data from single exposures are replaced by the data of adjacent exposures. A movie (m5.mpeg) 
showing the evolution of AIA 193 A is available online. 
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Peak Intensity (arbitrary unit) Velocity (km/s) Width (km/s) RB Asymmetry [70-130 km/s] 



2.48 3.01 3.55 4.09 4.63 5.16 5.70 -42 -28 -14 14 28 42 56 62 68 74 80 86 92 -0.21 -0.14 -0.07 0.00 0.07 0.14 0.21 
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Fig. 6. — Same as Figure[T]but for Fe XII 195. 12A in the 2007 Jun 5 observations. The rectangular region marks the 
locations where line profiles are averaged and presented in Figure [T3l 



Fe XII 195. 12A line since the exposure time used in 
the associated observations is too short so that only the 
strong Fe XII 195. 12A line has enough S/N to allow a 
reliable RB asymmetry analysis to individual profile. 
Since the blend Fe XII 195. 18A sits a t the red wing 
of Fe XII 195. 12A dYoung et alJl2009h . any blueward 
asymmetries detected by our RB p technique are not 
caused by this identified blend. 

The pre-emption parameters are presented in the 
first row of each figure. The scanned regions for all 
rasters are almost the same for almost every event. The 
only exception is the 2007 May 19 event shown in Fig- 
ure [21 where we can clearly see that the observed re- 
gion in the pre-emption phase is about 50" smaller in 
solar Y, compared to that in the eruption phase. 

3. Flows 

We found various types of flows in our observa- 
tions. In the following we mainly investigate proper- 
ties of three types of outflows associated with coronal 
dimmings and CME or EUV jet eruptions. 

3.1. High-speed outflows in dimming regions 

Coronal dimmings are clearly seen from the in- 
tensity images presented in Figures OB The dim- 



ming regions are characterized by a blueshift of 10- 
40 km s -1 , a notable phenomenon in the Hinode era 



dHarra et all 120071: 1 Jin et al.1 120091: lAttrill et al.1 l2010l: 
Chen et al.ll2oToaf lHarra et alJl201 lab . Enhancement 
of the line width in dimming regions has also been 



or me line widtn in dimming regions nas also been 
reported bvlMcIntoshl (|2009h . IChen etaD feOlOah and 
iDolla & Zhukovl d201 1 ) and it is very clear from Fig- 
ures HEl The significant blueshift and enhanced line 
width are similar to those found at the weak-emission 
boundaries of ARs (e.g., iMarsch etaD 120041 [20081; 



Harra et al.l 120081: IPel Zannal 120081: Ibel Zanna et al 



2011; Doschek et al. 2007, 2008; Tripathi et al. 2009; 



He et al. 2010a; Murrav et al. 2010; Brooks & Warren 
201 l l ilWarren et al.l201 llilBradshaw et al.l201 llifScott & Martensl 
201 ll: lYoung et al.ll2012l: &er et al.ll2012|; lHara et alJ " 



20081: iDe Pontieu et al. 1 [20091: iDe Pontieu & Mclntos 
201Q|: | McIntosh & De Pontieu 2009a b; Mc intosh et al 



201ll:IPeteii201Ql:lBrvans et al.l2010l:IUgarte-Urra & Warrenl 



201 ll:lMartinez-Svkora et al.l201 llfean et al.l201 lal^T 
FiguresHHalso reveal a significant blueward asym- 
metry in dimming regions. We note that the blue- 
ward asymmetries on maps of RB5 asymmetry for 
the 2006 Dec 14-15 observation s, which were pre- 
sented by Mci ntosh et al.l d2010h . are not so promi- 
nent as those in the RBp asymmetry maps in our 
Figure Q] This is because of the underestimation 
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Fig. 7.— RB asymmetry profiles (bottom) of the Fe x 184.54A, Fe xn 195.12A, Fe xm 202.04A and Fe xiv 274.20A 
line profiles (top) averaged over the square marked in Figure[T] Top: The observed spectra and measurement errors are 
shown as the diamonds and error bars, respectively. The green lines are single Gaussian fits. The two dashed red lines 
in each panel represent the two Gaussian components and the solid red line is the sum of the two components. The 
velocity (v) and exponential width (w) derived from the single (SGF) and double (I st 12 nd for the two components) 
Gaussian fits are shown in each panel. Also shown is the intensity ratio of the secondary component to the primary 
one (i2/il). Bottom: the black and blue lines represent RB profiles for RB p and RB^, respectively. Error bars indicate 
the errors propagated from the measurement errors. The peak relative intensity (i), velocity (v), and 1/e width (w) are 
shown in each panel. 
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Fig. 8. — First & second columns: RB asymmetry profiles (bottom) of Fe XIII 202. 04A line profiles (top) averaged 
over regions 1 & 2 marked in Figure 12 Third & fourth columns: RB asymmetry profiles (bottom) of the averaged 
Fe XII 195. 12A line profiles (top) in dimming regions at 02:51 and 03:18 as shown in Figure |U The line styles and 
denotations of parameters are the same as in Figure [7] 
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or me degree or asymmetry by trie K±$ $ tecnnique 
(Tianet al. 2011$. The presence of obvious blueward 
asymmetries suggests that the line profiles in dim- 
ming regions probably contain a highly blueshifted 
secondary component besides the primary compo- 
nent. Such a scenario is similar to that of the chromo- 



spheric network and AR edges (e.g..lHara et al. 200: 
De Pontieu et al.ll2009l: Be Pontie u & McIntost jEoT 
Mcintosh & De Pontieul2009 alibi : iMcIntosh et al.l201 1 
Peterj |2010t iBrvans et al.l 12010; iTian et al.l |2011a||c|; 



Ugarte-Urra & Warrenl 1201 ll: iMartinez-Svkora et al 
201 1|) . and thus we can perform a similar analysis of 
the line profiles. For further analysis, we only selected 
those locations where the average RB p asymmetry in 
the velocity range of 70-130 km s _1 is smaller than 
-0.03 (obvious blueward asymmetry) and the signal to 
noise ratio of the profile (defined as the ratio of the 
peak and background intensities) is larger than 8. The 
R B p -guided double Gaussian fit algorithm (see details 



Tian et al.ll2011cl) was then applied to the profiles at 



these locations. After the double Gaussian fit, we took 
the spectral position of the primary component as the 
line centroid and calculated the RB p> asymmetry pro- 
file. 

In Figure |7]&[U we present several examples of the 
observed and fitted line profiles and the correspond- 
ing RB p and RB d asymmetry profiles in dimming re- 
gions. By comparing the observed profiles with the 
different fitting profiles, we can clearly see the better 
performance of the double Gaussian fits and the devia- 
tions of the observed profiles from the single Gaussian 
fits. The RB asymmetry profile is basically the dif- 
ference between the emission of the two wings as a 
function of spectral distance (expressed in the veloc- 
ity unit) from the line center. Here a negative value 
means that the blue wing is enhanced with respect to 
the red wing at a certain spectral distance. The relative 
intensity and velocity of the secondary component can 
be derived from the peak of the RB asymmetry pro- 
file. And the 1/e width of the RB asymmetry profile is 
taken to approximate the width of the s econdary com- 
p qnent (see details inlTian et"aDl2Ql lcb . As discussed 



Tian et al.l (1201 lcb . the blend Si VII 274. 18A should 



in 

not have an important influence on the results of our 
RB asymmetry analysis and double Gaussian fit for 
Fe xiv 274. 20A since the two lines are very close to 
each other and the maps of the Fe xiv 274. 20 A line pa- 
rameters resemble those of the clean Fe xiv 264. 78A 
line. Moreover, as mentioned in the following, the 
contribution of Si VII 274. 18A to the total emission 



is at most 5.4%. The Fe xil 195. 12A line is blended 
with Fe XII 195. 18 A, which sits at the red wing of 
Fe xil 195. 12A so that any enhancement on the blue 
wing of the line profile is not caused by this identi- 
fied blend. Thus, the blueward asymmetries of the 
Fe xil 195. 12A line profiles we observed here are real. 
The Fe x 184.54A line is not as strong as the other 
three lines and there is a weak F e xi 184.41 A line at 
the blue side ( Brown et a l. 2008), making it difficult 
to derive the real degree of profile asymmetry. How- 
ever, the Fe xi 184.41 A line is about 210 km s _1 away 
from the Fe X 184.54A line and in normal conditions 
the two lines show up as two distinct peaks. So we 
believe that the enhancement between the two lines 
are due to the Fe X 184.54A emission from a high- 
speed (~100 km s _1 ) outflow. As can be seen from 
Figure [71 we usually found that the blueward asym- 
metry is present at different coronal temperatures and 
that the velocity of the high-speed upflow does not 
show a dramatic change with temperature. However, 
we are aware that an accurate velocity determination 
is beyond the ability of the EIS instrument because of 
the large instrumental width and the complication by 
blends. 

The ridge of enhanced line width for the 2007 May 
19 11:41-15:23 sca n, as shown in Fig ure [2l was pre- 
viously reported by IChen et al.l ( 2010a|) . They found 
that this ridge corresponded to the outer edge of the 
dimming region. In the Dopplergram, blue shift of 
~ 10 km s " 1 seems to be pres ent in the southern part of 
the ridge dChen et al.ll2010a|) . While the northern part 
shows net r ed shif t, which seemed to be omitted by 



Chen et al.l (1201 0a|) . From the map of RB asymmetry, 



we can see that this ridge is also charactered by clear 
blueward asymmetry. Typical line profiles in both the 
southern and northern parts of the ridge are presented 
in Figure [U Results of our RB asymmetry analysis in- 
dicate that the high-speed upflow is present along the 
whole ridge and that its speed might be ~ 100 km s _1 . 

We derived the parameters (relative intensity, ve- 
locity, exponential width) of the high-speed secondary 
component by using the three methods (double Gaus- 
sian fit, RBp and RB^) and in Figure we present 
the histograms of these three parameters derived from 
Fe XIII 202.04A line profiles in the 2006 Dec 15 01 : 15- 
03:30 observation (use profiles within the red contours 
shown in Figure [I]). We can see that the relative inten- 
sity is usually around 10% and can sometimes reach 
more than 30%. The velocity is usually in the range 
of 50-150 km s _1 and its distribution peaks around 
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Fig. 9. — Histograms of the relative intensity (A), velocity (B), and exponential width (C) of the secondary component, 
as derived from double Gaussian fit (black) and RB asymmetry analysis (red/blue for RB p/RB^) for Fe XIII 202. 04 A 
in the 2006 Dec 15 01:15-03:30 observation. Panel (D) shows the histograms of the xl values of the single (red) and 
double (black) Gaussian fits, as well as the ratio of the two (blue). Panels (E)&(F) show the relationship between 
the intensity/Doppler shift and the exponential width as derived from single Gaussian fit. Panels (G)&(H) present 
the relationship between the Doppler shift/exponential width derived from single Gaussian fit and the average RB p 
asymmetry in the velocity interval of 70-130 km s _1 . 
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90 km s _1 . The distribution of the exponential width 
peaks around 55 km s _1 , which is comparable to the 
width of the primary component. The distributions of 
the Xr f° r b° m the single and double Gaussian fits 
peak at values smaller than unity, which may result 
fro m the overestimation of the E IS measurement er- 
ror dPeterfeOldlTian et al.ll2011ch . However, from the 
Xr ratio between the double and single Gaussian fit we 
can see that the double Gaussian fit does better than the 
single Gaussian fit for these asymmetric line profiles. 

In Figure [9] we also show the relationship between 
the intensity/Doppler shift and the exponential width 
as derived from single Gaussian fit, and the rela- 
tionship between Doppler shift/exponential width de- 
rived from the single Gaussian fit (SGF) and the av- 
erage RB p asymmetry in the velocity interval of 70- 
130 km s _1 . There seems to be a weak anti-correlation 
between the SGF intensity and line width, which is 
consistent with previous result that the intensity and 
line width sho w negative corre lations in loop foot- 
point regions ( iScott & Martensl l201l|) although the 
correlation turns into positiv e when considering the 
whole AR dLi&Dingll2009h . Panel (F) reveals an 
obvious correlation between the SGF Doppler shift 
and line width. The calculated correlation coeffi- 
cient is -0.59 if using all data points inside the red 
contours. Similar correlation was also found at AR 



boundaries dDoschek et aDl2QQ7l [2008). As proposed 
by iDoschek et al. d2008|) . this correlation may sug- 
gest that the profile is composed of multiple compo- 
nents. Indeed, we find striking correlations in panels 
(G)&(H) of Figure with a correlation coefficient of 
0.67 in (G) and -0.58 in (H). Such correlations strongly 
suggest that the clear blue shift and enhanced line 
width in dimming regions are largely caused by the 
blueward asymmetries. T he growth of Al fven wave 
amplitude, as suggested by iMcIntoshl ( 2009]) . may be 
an additional reason for the enhancement of the line 
width. The fact may be that there is a faint high- 
speed upflow superimposed on a strong and almost 
stationary (or slightly shifted) background in the LOS 
direction. This scenario is also similar to the inho- 
mogeneiti es of flow veloc i ties al ong the LOS as pro- 
posed by iDolla & Zhukovl (|201 1|) and would naturally 
produce blueward asymmetric line profiles. A single 
Gaussian fit to the total emission line profile would 
yield a blueshift and enhanced line width, as com- 
pared to the line profile of the background emission. 
When the relative intensity of the high-speed upflow 
component becomes larger, the blueward asymmetry 



becomes more obvious and we will obtain a larger 
blue shift and line width if applying a single Gaussian 
fit. The intrinsic assumption of a single Gaussian fit 
is that everything is moving at the same bulk speed, 
which is obviously not the case in coronal dimming 
regions. Thus, our analysis implies that previous re- 
sults based on a single Gaussian fit can not reflect the 
real physical processes and thus may need to be recon- 
sidered. First, the outflow speed i n dimming regions 



is perhaps not around 20 km s MHarra et al] 12007 
iTF 



rtaj] 

Jin etalJ 120091: lAttrill et alJ 1201(1 IChen et al .112010a' : 



Harra et al.ll2011a|) . but can easily reach ~100 km s _1 
in the lower corona. Second, the enhanced line width 
is not pure ly due to the in crease of the Alfven wave 
amplitude (iMcIntoshll2009|) . but is largely contributed 
by the superposition of different emission components. 

We note that the properties of these high-speed 
outflows are very si milar to those we found previ- 
ously in AR edges dTian et al.1 1201 lc|) . This sim- 
ilarity suggests that the outflows in both regions 
may result from a similar process, e.g., h eating 
in the lower atmosphere dl)e Pontieu et al. 2009 : 



Mcintosh & De Pontieul 12009 bt lHansteen et al] 12010: 



Song & Vasyliunasl l2011[) . Magnetohydrodynamic 
simulations have shown that magnetic reconnection 
is an efficient mechanism to produce high-speed 
outflows (jets) in the lower solar at mosphere (e.g., 
Ding et al.1 Hoill: iRoussev et all 1200 lb . As the mag- 
netic field lines opened up by CMEs, rapid multi- 
thermal upflows produced by both the pre-existing and 
CME-induced impulsive heating at the lower part of 
the erupted loops are guided by the field lines into 
the transiently opened corona. These outflows may 
serve as an important source of materials to refill the 
corona. From Figures [TJ&14] we can see that the blue- 
ward asymmetries were strongest within a few hours 
after the flare peak time, indicating that the high-speed 
outflows were strongest right after the erupted materi- 
als left the Sun. As the dimmings gradually recovered 
and the magnetic fields began to close down again, the 
outflows became weaker. Such a resul t is consistent 
with the finding of iMiklenic et al.l d201l|) that the mass 
loss occurs mainly during the period of strongest CME 
acceleration. 

Through joint ima ging and spectroscopic observa - 
tions of the corona. iMcIntosh & De Pontieul ( 2009a) 



De Pontieu & Mcintosh! d2010h and lTian et allBoiTalcT) 
have suggested that the secondary emission compo- 
nent found at AR edges is caused by high-speed repet- 
itive upflows in the form of upward propagating distur- 
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bances (PDs) in EUV and X-ray imaging observations. 
Similarly, we think that the highly blueshifted compo- 
nent found in spectra of dimming regions should ex- 
hibit as PDs in imaging observations. The SDO/AIA 
observations, with a high S/N (especially in the 171 A 
& 193 A passbands) and high cadence, might be able to 
reveal such weak PDs. The 201 1 June 21 observations 
were done in the SDO era and from the associated 
movie (m4.mpeg) of Figure |4] we can clearly see the 
evolution of the dimming boundary, which is likely to 
be assoc iated with the success ive disappearance of the 
moss (Mcintosh et al. 2QQ^). And, there seems to be 
weak outward PDs along legs of the opened coronal 
loops. However, both the LOS effect and the signifi- 
cantly reduced emission associated with these opened 
field lines make it difficult to study the PDs quantita- 
tively. Part of these outflows may experience further 
acceleration at higher layers, overcome the gravity, 
and eventually become the solar wind stream along 
the transiently opened field lines, which may serve 
as an additional momentum source for the associated 
CME dMdntosh et alJl2010h . We also noticed that the 
ascending post-flare loops revealed by the movie were 
clearly observed by EIS (last row of Figure |4]), show- 
ing a blue shift of ~5 km s _1 . 

3.2. Temperature-dependent outflows 

A temperature-d ependent outflow was reported by 



Imada et al.l (I2007|) . who found that the flow speed 



increases from ~10 km s _1 at log (77K)=4.9 to 
~150 km s _1 at log (77K)=6.3 in a dimming region. 
However, we found that this temperature-dependent 
outflow is not in but immediately outside the deep- 
est (darkest in the intensity image) dimming region. 
More interestingly, we found that our RB p technique 
can identify this temperature-dependent outflow. This 
event was located around (x=135 // ,y=-35 // ) in Figure [3] 
and we can see that it is associated with a small patch 
of redward asymmetry. This is easy to understand 
since the outflow component is much stronger than the 
background emission component in Fe XIII 202.04 A, 
as can be seen from Figure [3f of llmada et al.l d2007|) . 

We further identified several temperature-dependent 
outflows in the 2006 Dec 13 01:12-05:41 observation. 
These temperature-dependent outflows are associated 
with the small patch of redward asymmetry around 
the location of (x=130",y=-85"), (x=175",y=-95"), 
(x=375",y=-135") and (x=135",y=-35"), respectively. 
All of these temperature-dependent outflows are not in 
but immediately outside the deepest dimming region. 



As an example, Figure [TO] shows the line profiles of the 
temperature-dependent outflow around (x=175 // ,y=- 
95"), region 1 in Figure [3] A single Gaussian fit 
seems to be adequate to derive the outflow veloci- 
ties for the emission lines formed at a temperature of 
log (77K)<6.0. For emission lines formed at higher 
temperatures, we see clear indications of two well- 
separated components in the line profiles. Thus, we 
applied a double Gaussian fit to these line profiles and 
the Doppler shift of the highly blueshifted compo- 
nent (denoted as v2 in Figure [TO]) should represent 
the outflow velocity at the corresponding tempera- 
ture. The small-velocity component, whose velocity 
is denoted as vl, is likely to be the nearly stationary 
background emission of the corresponding ion. Since 
we are mainly interested in the velocity of the highly 
blueshifted component and this component is often 
stronger than the background component, the outflow 
velocity derived from the double Gaussian fit should 
be highly reliable. We can see that the temperature 
variation of the outfl ow velocity shown i n FigurefTOlis 
similar to Figure l6l of llmada et al.l (|2007|) . 

We also identified temperature-dependent outflows 
in the 2011 Jun 21 observation. Unfortunately, the 
exposure time used in this observation is only 9 sec- 
onds so that only the Fe xn 195. 12A line can be 
used for asymmetry analysis without any temporal 
or spatial binning. From Figure H] we can see en- 
hanced redward asymmetries surrounding the region 
of significant blueward asymmetry (dimming region). 
The blend of Fe XII 195. 18A certainly contributes 
to these redward asymmetries. However, if we spa- 
tially bin line profiles of various lines at several ad- 
jacent pixels there, we start to see clear signatures of 
temperature-dependent outflows. Thus, both the blend 
of Fe XII 195. 18A and the temperature-dependent out- 
flows are causing these redward asymmetries. As an 
example, we present in Figure[TT]the line profiles aver- 
aged over the small rectangular region marked in Fig- 
ure HI Similar to Figure [lOl single or double Gaussian 
fits are applied to these line profiles to obtain the out- 
flow velocities at different temperatures. 

From Figures [3jSj4l we can also see the enhancement 
of the line width at locations where the temperature- 
dependent outflows are found. These enhanced line 
widths, as derived from single Gaussian fits, are actu- 
ally caused at least partly by the superposition of the 
relatively weak background emission component and 
the strong outflow component for Fe xin 202. 04 A and 
Fe XII 195. 12A, as demonstrated in FigurefTOland Fig- 
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Fig. 11. — Same as Figure[K)]but for the rectangular region marked in Figure[4] 
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ure[TTJ respectively. 

Note that there are several blends of the He II 256. 32 A 
line. However, the He II 25 6. 32 A line usually domi- 
nates and contributes more th an 80% of the tota l emis- 
sion in disk observations (Youn g et al.l 120071) . The 
Fe VIII 185.21 A is blended with Ni XVI 185. 23 A but 
the blend should not have a large impact on the de- 
rived velocity, since the two lines are very close to 
each other and the latte r is much weaker than the for- 
mer (I Young et al.ll2007|) . The enhancement of the blue 
wing in the Fe XV 284. 16A line profile, at around 
283. 95 A in Figure fTTl seems to be caused by the weak 
blend Al ix 284.03 A( Young et alJ [2007b and should 
not impact the derived velocity of the very strong out- 
flow component significantly. 

The fact that these temperature-dependent out- 
flows are found outside the (deepest) dimming re- 
gions suggests that these outflows are different from 
the high-speed outflows we described in the previ- 
ous section. The temperature-dependent nature of 
these outflows resembles that of gentle (as opposed 
to explosive) chromospheric evaporation flows. Gen- 
tle chromospheric evaporation can be driven by low- 
flux (<10 10 ergs cm -2 s _1 ) nonth ermal electron 



beams in the flare impulsive phase (e.g jMiUigan et al 
l2QQ6b b or thermal conduction in the flare decay phase 



re.g jAntiochos et al.|[l978l: iBerlicki et al.ll2005l) . The 
temperature-dependent outflows we present here are 
away from the flare sites so that they may not be di- 
rectly related to the associated flares at first thought. 
However, we can not exclude the possibility that some 
magnetic field lines there are connected to the flare 
sites and that nonthermal electron beams or enhanced 
thermal conduction resulting from the flares cause the 
evaporation flows. For the 2006 Dec 12-13 event we 
may exclude the possibility of nonthermal electron 
beams since most temperature-dependent outflows 
were identified after the flare peak time. It is also pos- 
sible that interactions between the opened field lines 
in the dimming region and the closed loops in the sur- 
rounding plage region produce low- flux nonthermal 
electrons or/and enhanced thermal conduction which 
will then generate the evaporation flows. 

The difference between the high-speed outflows in 
dimming regions and the temperature-dependent out- 
flows immediately outside the (deepest) dimming re- 
gions can also be understood in the sense of driving 
force. The former are perhaps driven by magnetic re- 
connection in the chromosphere or TR. These multi- 
thermal outflows usually do not show obvious temper- 



ature dependence since the acceleration by magnetic 
force finishes at locations very close to the reconnec- 
tion site (e.g., lYokoyama & Shibatalll995h . As dis- 
cussed above, the latter seem to be evaporation flows 
which are driven by pressure gradient force. Thus, 
the acceleration wil l continue as long as a pressure 
gradient exists (c-g -, KamioetaL|[2009; Shi mojo et al 



200 ll: lJudge eTSfeoi2h . The dramatic increase of the 
flow speed from log (77K)<6.0 to log (77K)>6.0, as 
shown in Figu re ITOl is s imilar to that of the event 
analyzed by llmada et al.1 (12007b and may be caused 
by a st eep pressure (temp erature) gradient at a certain 
height (llmada et alJl201lh . 



3.3. Highly blueshifted component representing 
the ejecta emission 

No pronounced dimming was recorded by EIS for 
the 2011 Feb 14 and 2007 Jun 5 observations. How- 
ever, clear plasma ejections (CME or EUV jet) were 
observed and the associated line profiles clearly ex- 
hibit two well- separated components. The eruption of 
the CME loop on 201 1 Feb 14 was clearly revealed in 
the AIA 304A, 17lA & 193Apassbands (see m5.mpeg 
showing the evolution of the 193 A passband, the green 
box indicates the FOV of EIS observation). For imag- 
ing observations of the 2007 Jun 5 jet, we refer to 
YangetalJ d201lh . Fi gures fT2bfT3l show an example 
of line profiles associated with the CME ejecta and the 
EUV jet, respectively. It is clear that emission lines 
formed at coronal and TR temperatures clearly exhibit 
two well- separated components, an almost stationary 
component accounting for the background emission 
and a highly blueshifted component representing emis- 
sion from the erupting material. The Doppler veloci- 
ties of the two components can be easily calculated 
through a double Gaussian fit. 

From Figure[l2]we can see that the highly blueshifted 
component has a velocity of ~220 km s _1 and the ve- 
locity does not change significantly with temperature. 
This velocity should represent the LOS velocity of the 
CME ejecta (expanded loops) at 19:29 on 2011 Feb 
14. The POS component of the ejecta velocity was 
estimated to be ^200 km s _1 from simultaneous AIA 
193 A observations. Combining the two velocity com- 
ponents, we can calculate the real speed of the ejecta 
at 19:29, which is ~300 km s _1 . Unfortunately, the 
fast EIS scans only focused on the same region so that 
the oblique propagation of the CME ejecta was not 
fully tracked. But in principle one should be able to 
track the complete evolution (not only the POS com- 
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Fig. 12. — Line profiles averaged in locations marked in Figure [5] The diamonds, error bars, line styles and colors are 
the same as in Figure [7] The peak intensities (il & i2), velocities (vl & v2) and exponential widths (wl & w2) of the 
two Gaussian components are shown in each panel. 



ponent) of CMEs by using simultaneous imaging and 
spectroscopic observations. We noticed that clear line 
splittings interpreted as filament or plasmoid erup- 
tions at similar s peeds wer e previous l y repo rted by 
Harra & Sterlmgl (12003b and lAsai et all d2008h . 

From Figure [13] we can see that the blueshifted 
component, which is apparently associated with the 
EUV jet, has a speed of ~170 km s _1 . The fact that 
we do not see an obvious temperature dependence of 
the flow speed suggests that this EUV je t is produced 
by re connection instead of evaporation ( Kamio etaD 
2009h . Since the PQS compon ent of the jet speed is 
~145 km s _1 dYang et al.ll201lh . the real speed is cal- 
culated to be ~223 km s _1 . From Figure [13] we can 
also see that the intensity ratio of the nearly station- 
ary component and the highly blueshifted component 
(il/i2) increases from 0.90 at log (77K)=4.7 to 5.85 at 
log (77K)=6.2. Such an increase should be directly re- 
lated to the difference in the temperature distribution 
(DEM, discussed below) of the background emission 
and the jet emission. It is likely that the decrease of 
blue s hift with temperature, as derived by I Yang et aD 
(1201 ll) from a single Gaussian fit, is in fact caused by 
this increase of intensity ratio and thus can not reflect 
the real physical process. In addition, the very large 



no n-thermal velociti es (~ 100-400 km s -1 ) reported 
by iKim et al.1 d2007h and lYang et al.1 d201lh through 
single Gaussian fits are also likely to be caused by the 
effect of the superposition of the two (background and 
jet) emission components. 

Another feature that is worth noting is the inverted 
Y- shape structure at one footpoint of the erupted loop 
(the fifth row of Figure [5]) and the base of the jet 
(the second row of Figure @). A net red shift and 
enhanced line width are found at the base of the in- 
verted Y- shape structure in each case. Red shifts have 
been previously reported at the ba se o f a polar jet 
and an AR jet bv lKamio etaD d2007h and lChifor etaT 



(2008a), respectively. It is likely that they are caused 
by the downward propagating reconnection outflows 
which collide with and comp ress the underlying loops 
dYokoyama & Shibatalll995|) . And the enhanced line 
widths perhaps result at least partly from the flow in- 
homogeneities in this process. It seems that significant 
redward asymmetries are also found at the bases of the 
inverted Y-shape structures. But they are complicated 
by the blend Fe XII 195. 18 A. However, by averaging 
profiles of some weaker but clean lines over several 
spatial pixels at the bases of the inverted Y-shape struc- 
tures, we do see signatures of redward asymmetries 
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Fig. 13. — Same as Figure[l2]but for the rectangular region marked in Figure[6] 



(not shown here) which may be caused by the down- 
ward flows. 

At the end of this section, we would like to discuss 
various types of flows found in the 2006 Dec 12-13 
observations . For detail e d desc r iptions of thi s even t, 



we refer to lKubo etal.1 d2007l). Izhangetal 



AsaietalJ d2008l) . Ijing et :alJ I d2008l) and Em 



( 2007 ). 
( 2011 ). 



From Figure [3] we can see that before eruption (19:07- 
23:46) the loop footpoint regions are characterized by 
clear blue shift, enhanced line width and obvious blue- 
ward asymmetry. Such results indicate the presence of 
a weak high-speed upflow sup erimposed on the nearly 
stationary background ( e.g., iDe Pontieu et alJ 20091: 
Mcintosh & De Pontieul \l0O9sk iTian et al.1 l2011allcl : 
Martinez-Sykora etaD 201 1|) . During eruption (01:12- 
05:41) there was a clear expansion of the regions with 
clear blue shift, enhanced line width and obvious blue- 
ward asymmetry. Line profiles in the newly formed 
dimming regions (e.g., profile of Region 3 shown in 
Figure [14]) are obviously blueward asymmetric and 
thus are similar to those of loop footpoint regions, 
indicating the presence of rapid upflows along the 
opened field lines as discussed in Section 3.1. The 
elongated blueward-asymmetry feature was observed 
in the flare impulsive phase and it was not located in 
the dimming region. The line profiles there (e.g., pro- 
file of Region 2 shown in Figure [Pfl) clearly reveal two 



components and the high-speed component is likely to 
be associated with the initial removal of the magnetic 
loop. From Figure[3]we can also see patches of signifi- 
cant redward asymmetries around the flare site (around 
x=330 // ,y=-80 // ). FigurefT4lshows profiles of two lines 
in a small region (Region 4 marked in Figure [3]) and 
the redward asymmetries are clearly revealed in both 
the single Gaussian fits and RB p asymmetry profiles. 
Such asymmetries clearly indicate that the enhanced 
nonthermal broadening i s caused by the superposition 
of flows, i.e., turbulence (Milligan 201 lh. We found a 
net red shift of ^30 km s _1 for almost all strong lines 
used in this observation. This multi-thermal down- 
ward motion is perhaps driven by both the cooling 
of the flare plasma and the overpressure of the flare 
plasma relative to the underlying atmosphere. 

4. Plasma diagnostics 
4.1. Dimmings 

Using CDS data lHarrison & Lvonsl d2000h and 
Harrison et al.1 ( 2003|) made first efforts to diagnose 
the electron densities and calculate the mass losses 
of dimming regions observed at the limb. They con- 
cluded that the reduced emission in dimming regions is 
an effect of mass loss rather than temperature change. 
They also mentioned the importance of calculating 
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Fig. 14. — Similar to Figure [7] but for regions 2-4 marked in Figure [3j For the averaged profiles in region 4 only the 
single Gaussian fits and the RB p asymmetry profiles are plotted. 



mass losses of on-disk dimming regions in the context 
of space weath er forecas t . Base d on static solar atmo- 
sphere models, I Jin et al.1 (12009b also tried to calculate 
the mass losses of two dimming regions by using EIS 
observations. 

Here we extend these previous investigations and 
make efforts to diagnose the electron density, temper- 
ature and mass loss for the dimming regions we study 
in this paper. By comparing the intensity images of 
various emission lines before (19:20-21:35) and af- 
ter (01:15-03:30) the eruption for the 2006 Dec 14- 
15 event (Figure [15]), we can clearly see the occur- 
rence of dimming at all temperatures. We then aver- 
aged line profiles over the regions where the intensity 
was reduced more than 20% and calculated intensities 
of different lines before and after the erupti on. By us 



ing th e routine chiantijdem.pro (also used by lLee et al 



201 lh available in SolarSoft and assuming a constant 



pressure of 10 lb cm -3 K, we obtained the differential 
emission measure (DEM) curves at the pre-emption 
phase and of the dimming region. Here a double Gaus- 
sian fit was applied to the Fe xi 188. 23 A line profiles 
to derive the input line intensities since Fe XI 188.23A 
is partly blended w ith the strong Fe XI 188.30A line 
(You ng et aDE0 07). While we applied single Gaussian 
fits to the profiles of other selected lines. This means 



that the blends to Fe xn 195. 12A and Fe xiv 274.20A 
were simply neglected since their contribution to the 
total emission is of the order of 5% or less in outer 
parts (low density) of ARs (see below and lYoung et al.l 
l2009h . The influence of the blend Al ix 284.03A 
is also negligible since it is at the far wing of the 
Fe XV 284. 16A line and is very weak in AR condi- 
tions. The He II 25 6. 32 A line is blended with sev- 
eral other higher- temperature lines (Fe XIII 25 6. 42 A, 
Fe XII 256.41A, Si x 256.37A). We included this line 
for our DEM analysis since it is the only strong EIS 
line formed in the lower TR and it contributes more 
than 80% of the t otal emission in disk observations 
d Young et alJ 120071) . Moreover, the blends are all sit- 
ting at the red wing of the He II 256. 32 A line profile 
and their spectral distances from the He II 256. 32 A 
line center are ~60-120 km s _1 . We note that in such 
a case our single Gaussian fit algorithm mainly fits the 
core and blue wing of the line profile, which is pri- 
marily the emission of He II 256. 32 A rather than the 
blends. Note that the weak high-speed outflow was 
not considered in the DEM and the following density 
diagnostics since it only contributes a few percent to 
the integrated intensity of the average line profiles. We 
can see from Figure [16| that the main difference is the 
reduced emission at high temperatures (log (77K)=6.1- 
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Fig. 15. — Intensity maps of different emission lines for the two scans starting at 19:20 on 2006 Dec 14 and 01:15 on 
2006 Dec 15. The red contours outline regions where the intensity was significantly reduced (^20%). 
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Fig. 16. — DEM curves for the two scans starting at 
19:20 on 2006 Dec 14 (upper panel, before eruption) 
and 01:15 on 2006 Dec 15 (lower panel, dimming). 
The DEM curve before the eruption is overplotted as 
the dashed line in the lower panel. 



6.3) in the dimming region. This result seems to sug- 
gest that in this event a significant portion of the cool 
TR materials did not escape when the magnetic field 
lin es opened up. Simila r result has also been obtained 
by iRobbrecht & Wangl d2010h based on EUV imag- 
ing observations. We have to mention that the lower- 
temperature part of the DEM curves is less constrained 
due to the lack of many cool lines in the EIS observa- 
tion. Observations of IRIS, which will be launched in 
2012, are thus crucial since several strong cool lines 
are included in its spectrum. 

There are several density sensitive line pairs used 
in the first (15:11-16:01) and last scans (10:29-11:19) 
of the 2006 Dec 14-15 observation. Thus, we can 
compare the densities of the dimming region with pre- 
emption densities. Line pairs Fe xn 186.88A& 195. 12A 
and Fe xm 203. 82A & 202.04A were chosen for the 
density diagnostics. Fe XII 186.88 is actually a self 
blends of two Fe xn lines (186.85A, 186.88A) and the 
Fe XII 195. 12A lin e is blended with Fe xn 195. 18A 
dYoung etalJl2007h . The S xi 186.84A line usually 
contributes no more than 5% of the Fe xil 186.88A 
feature and thus w as ignored in our calculation 
dChifor etalJl2008ah . Fe xm 203. 82A line is a self 
blends of two Fe xm lines (203.79A, 203. 82 A). We 
calculated the line ratios and the associated standard 
deviations in the obvious dimming region (defined as 
the region with an intensity reduction larger than 20%) 
of the last scan and in the corresponding region of the 
first scan. The theoretical relationships between the 
line ratios and densities, as extracted from CHIANTI 
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Fig. 17. — Relationship between electron density and 
line ratio, as derived from the CHIANTI database. The 
measured values before the eruption (15:11 on 2006 
Dec 14) and during dimming (10:29 on 2006 Dec 15) 
are indicated by the solid diamonds and squares re- 
spectively. The vertical and horizontal bars indicate 
the standard deviations of the intensity ratios and un- 
certainties of the calculated densities, respectively. 



database (Dere et al. 1997; LandietaL 2006), are pre- 
sented in Figure [T7] The measured values before the 
eruption and during dimming are indicated by the solid 
diamonds and squares respectively. We can see that 
the average density changes from log (A^ e /cm _3 )=8.89 
to log (Af e /cm- 3 )=8.67 at log (77K)=6.1, and from 
log (Af e /cm- 3 )=8.70 to log (Aycm- 3 )=8.58 at log 
(77K)=6.2. This density decrease, together with the 
fact that the dimming is seen at all temperatures, 
strongly suggests that the dimming is an effect of den- 
sity decrease rather than temperature change. 

Following Harrison & Lvons (2000) and Uin et al 



(2009), we have attempted to estimate the mass losses 
in several well-observed dimming regions by using 
two different methods. The first method is just to 
multiply the density change and the emission vol- 
ume, which is simi lar to the Si X method used by 
lHarrison et al.1 d2003l) . The calculation process can be 
expressed as following: 



M = dNSLrriv 



(2) 



where M, 8 TV, S, L and m p represent the to- 
tal mass loss, change of the number density, area 
of the dimming region, depth of the dimming re- 
gion and proton mass, respectively. The density 
change SN can be calculated from the line pair 
Fe XII 186.88A & 195.12A, or estimated from the 
intensity change of Fe XII 195. 12A if the other line 
was not used in the scan of dimming. The use of 
the Fe xn 195. 12A line seems to be reasonable since 
from Figure [16] we can see that the most significant 
decrease as well as the DEM peak occur around log 
(77K)=6.1. The dimming area S is defined as the total 
area where the Fe xn 195. 12A intensity drops more 
than 20%, multiplied by the ratio of the total dimming 
area in simultaneous full-disk coronal images (EIT or 
AIA observations) and the dimming area observed by 
EIS. Assuming that the emission volume is as deep 
as it is wide, the depth of the dimmin g region can be 



calcul ated as y/S. As mentioned by lHarrison et al 



(2003), the mass calculated from this method should 
be considered to be a reasonable figure for a compar- 
ison with the mass of the associated CME, although 
the Fe xn 186.88A & 195.12A line pair can only de- 
tect density changes at a temperature of around log 
(77K)=6.1. 

The second method is similar to that used by 



Jin et al. (2009). We take the e mission heights of TR 
lines from the VAL3C model ( Verna zza et al.l Il981|) 
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and coronal lines from iMariska & Withbrod (I1978|) . 
then calculate the densities at these heights from an 
empirical relationship b etween height and density 
dCoxll200Ql:ljin et al.ll2QQ9h . Any density sensitive line 
pairs available in EIS observations are then used to de- 
rived the densities, which are then compared with and 
scale the model densities. Density changes at different 
heights can then be derived from intensity changes of 
corresponding EIS lines. The total mass loss is then 
expressed by 



M = y £SN(hi)S(h i )Sh i m p 



(3) 



where SN(hi) and S(hi) are the density change and 
dimming area respectively at the emission height of the 
/th line. 

The mass losses calculated from the two methods 
for several well-observed dimming events are listed 
in Table [2 It can be seen that the mass losses es- 
timated from different methods are 20% -60% of the 
CME mass cal culated from LASCO white light data 
djin et al.f2 009). Such results indicate that a significant 
part of the CME mass originates from the dimming re- 
gion. Thus, in principle the values of mass loss esti- 
mated from spectroscopic observations can be used to 
guide our identifications of CME/ICME sources. We 
have to mention that no strong TR lines were used in 
the 2011 Jun 21 observation so that the mass loss de- 
rived from the second method is likely to be underesti- 
mated. We think that the difference of the DEM curves 
before and during the dimming, as shown in FigurefT6l 
can in principle be used to derive the mass loss of the 
dimming region. However, it seems that we still have 
to make several assumptions for the dimming depth as 
well as the changes of density and temperature gradi- 
ent, which are hard to evaluate. Thus, we do not make 
an effort in this direction and leave it open for future 
investigations. 

Jin etal.1 12009) ade an effort to calculate the mass 



flux for the outflows in dimming regions and con- 
cluded that the total ejected mass is about one order 
of magnitude larger than the CME mass. The flow ve- 
locities they used were derived from a single Gaussian 
fit to line profiles in dimming regions and they are of 
the order of ~20 km s _1 . The density values were 
taken from static solar atmosphere models. If we as- 
sume that the mass refilling the corona comes from the 
high-speed outflows, we should expect to see an equiv- 
alence of the total mass supplied by the outflow and the 
mass loss in the corresponding dimming region. If we 



simply estimate the total mass supplied by the outflow 
as the product of mass flux density, area (S) and du- 
ration (t) of the dimming, we can have the following 
relationship: 



6NSL = nvSt, 



(4) 



From EIT or AIA observations, we roughly es- 
timated the duration of significant dimming as 14, 
14 and 10 hours for the 2006 Dec 14, 2006 Dec 13 
and 2011 Jun 21 events, respectively. If we take a 
value of 100 km s _1 for the speed (v), the density 
(n) of the outflow can then be calculated from Equa- 
tion |U The calculated densities of the outflows are 
listed in Table [2 These values are about two orders 
of magnitude smaller than the pre-emption densities 
at log (77K)=6.1 and should only be regarded as the 
lower limits since some of the outflows may over- 
come the gravity and become part of the solar wind. 
The lower limit of the mass flux density associated 
with these high-speed outflows is thus estimated to be 
about 1.67 x 10 -10 g cm -2 s _1 if using a density of log 
(Af e /cm- 3 )=7.0. While the average outflow mass flux 
density of the dimming region in the 01:15-03:30 scan 
on 2006 Dec 14 is estimated to be about l.OxlO -9 
g cm -2 s _1 if using results of single Gaussian fit, 
i.e., an average velocity of 12.5 km s _1 and a density 
of log (Af e /cm _3 )=8.67. These values are within one 
order of magnitude of the mass flu x density of type-I 



spicu les (^1.67 x 10~ 9 g cm 
1978b, tvoe-II spicule s (-1.5x10 
De Pontieu et alJE 



Pneuman & Koppl 



g cm s \ 
2011b. coronal rains (~1.14 x 10~ 9 



g cm 2 s \ Antolin & Rouppe van der Voortl 12012 ) 
and outflows in the quiet-Su n network (~1.6xl0 -9 
Tian et alJE009b , 



g cm 

We have to point out that the calculations of the 
mass flux of the high-speed outflows are based on 
the assumption that there are only two emission com- 
ponents and that the primary component is at rest. 
However, in fact we can not tell whether the primary 
component is really at rest or moving upward with 
a small velocity (e.g., ~10 km s _1 ) since EIS does 
not allow an absolute wavelength calibration. More- 
over, due to the large instrumental width we can not 
rule out the possibility of more than two components 
with each s lightly Doppler- shi fted with respect to each 
other (e.g JDoschek et alJEooik . 
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4.2. Erupted CME loop and EUV jet 
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Fig. 18. — DEM curves of the background (upper 
panel, stationary component) and ejecta (lower panel, 
highly blue- shifted component) for the rectangular re- 
gion marked in Figured The DEM curve of the back- 
ground is overplotted as the dashed line in the lower 
panel. 



In the case of erupted CME loops or EUV jets, since 
we can often unambiguously separate the ejecta emis- 
sion component from the background emission com- 
ponent, e.g., Figures [T2lSfT3l in principle we should be 
able to diagnose the plasma properties separately for 
each of the two components. Using the fitting results 
shown in Figure \13\ we calculated the DEM curves 
of the two components. From Figure [18] we can see 
that the ejecta (EUV jet) has more emission around 
log (77K)=5.5, as compared to the background. We 
can also see that the emitting materials of the EUV 
jet are almost equally distributed over the temperature 
range of log (77K)=5 .4-6.1. We have to mention that 
there was only one line with a formation temperature 
of log (77K)<5.5 (He II 256.32A) in our EIS observa- 
tions so that the low-temperature part of the DEM is 
not well constrained. Future joint observations of EIS 
and IRIS are thus highly desired. Unfortunately, there 
are no very cool lines (with a formation temperature 
comparable to that of He II 256.32A) in the 201 1 Feb 
14 observation so that we could not perform a reliable 
DEM analysis (with a wide temperature coverage) for 
the CME ejecta in this observation. 

The density sensitive line pair Fe XIV 264.78 A & 274.20A 
were included in both the 2007 Jun 5 and 2011 
Feb 14 observations. The Fe XIV 274.20A line 
is blended with Si VII 274. 18 A. The intensity of 
Si VII 274. 18A can be estimated from Si yn 275.35A 
since the ratio of the two is at most 0.25 dYoung et alJ 
l2QQ7h . We used another observation which included 
both the Si VII 275. 35A and Fe xiv 274.20A lines 



(ITian et al.ll2011q) and estimated that the contribution 
of Si VII 274.18A to the total emission is at most 5.4%. 
Thus, we simply ignored this blend. The theoretical 
relationship between the line ratio and density, as well 
as the measured values for the two components, are 
shown in Figure [I9j We can see that the measured 
values are in the density sensitive part of the theoret- 
ical relationship. The ejecta, which is an EUV jet in 
the 2007 Jun 5 observation and an erupted loop in the 
2011 Feb 14 observation, seems to have a larger den- 
sity compared to the background emission. The mea- 
sured densities are log (N e /cm _3 )=9.89 for the EUV 
jet and log (# e /cm _3 )=10.01 for the erupted CME 
loop. However, the uncertainties of the line ratios for 
the ejecta are very large. This is largely due to the fact 
that the ejecta component in the Fe xiv 274. 20 A line 
profile was too close to the edge of the spectral win- 
dow so that it was only partly resolved. The measured 
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Fig. 19. — Relationship between electron density and 
line ratio, as derived from the CHIANTI database. The 
measured values for the emission of the background 
and ejecta are indicated by the solid diamonds and 
squares respectively. The error bars indicate uncer- 
tainties of the ratios and densities as calculated from 
the 1-cr uncertainties of the double Gaussian fit inten- 
sities. Upper: the 2011 Feb 14 CME; Lower: the 2007 
Jun 5 jet. 



densities of the background components are consistent 
with (in the jet case) or about two times larger than 
(in the CME cas e) the normal AR d ensities at simi- 
lar temperatures dTripathi et al.ll2QQ8|) . To the best of 
our knowledge, this is the first time that the electron 
densities of EUV jets and erupted CME loops are mea- 
sured through an unambiguous decomposition of line 
profiles. Unfortunately, no density sensitive line pairs 
formed at lower temperatures were included in these 
two observations. 

For the EUV jet, the total mass of the ejected ma- 
terial can be calculated by multiplying the mass flux 
density, cross section area and life time of the jet. 
Taking a density of log ( J /V e /cm _3 )=9.89, a speed of 
223 km s -1 , a jet width o f 8" and a life time of 11 
minutes dYang et al.1 l201ll) . the total ejected mass is 
estimated as 5.0xl0 13 g. This is about two orders of 
magnitude lower than the typical value of CME mass. 

The mass of the erupted loop in the 201 1 Feb 14 ob- 
servation can be calculated by taking a density value of 
log (Af e /cm _3 )=10.01. After estimating the loop cross 
section area and loop length from the AIA 193 A im- 
age at 19:29, the mass of the erupted loop was esti- 
mated to be 2.5 xlO 14 g. Such a value is comparable 
to the lower limit of the typical CME mass. However, 
we have to bear in mind that the density value we used 
only represents the density of the emitting materials 
with a temperature around log (77K)=6.25. 

We have to mention that both the DEM and density 
calculations are based on ionization equilibrium. In 
the case of flows, this equilibrium might be destroyed 
(IPeter et al.ll2006|) . To quantify how the flows impact 
the results of temperature and density diagnostics, fur- 
ther numerical simulations are needed. 

5. Conclusion 

We have analyzed several data sets obtained by EIS 
during solar eruptions such as CMEs, coronal dim- 
ming s and EUV jets. We have mainly identified three 
types of flows and investigated the properties of them. 
We have also performed density diagnostics and DEM 
analyses for coronal dimmings, erupted CME loops 
and EUV jets. Our analyses suggest that spectroscopic 
observations can provide valuable information on the 
LOS kinematics and plasma properties of CMEs and 
EUV jets. 

Previous analyses based on single Gaussian fits re- 
veal significant blueshift and enhanced line width in 
the CME-induced dimming regions. However, our de- 
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tailed RB asymmetry analyses and RB -guided double 
Gaussian fits of the coronal line profiles clearly show 
blueward asymmetries in dimming regions, suggest- 
ing perhaps the presence of a relatively weak (~10% 
of the total emission) high-speed (~100 km s _1 ) up- 
flow component superimposed on a strong background 
emission component. This upflow component may re- 
sult from the impulsive heating in the lower solar at- 
mosphere. We have found that both the blue shift and 
line width correlate very well with the blueward asym- 
metry, suggesting that the significant blueshift and en- 
hanced line width are actually largely caused by the su- 
perposition of the two components. This finding sug- 
gests that a small portion of the plasma in the dimming 
region is flowing outward at a velocity of the order of 
100 km s _1 . Our plasma diagnostics of the dimming 
region suggest that dimming is mainly an effect of 
density decrease rather than temperature change. The 
mass losses in dimming regions have been estimated 
from two different methods and they are 20% -60% of 
the masses of the associated CMEs, suggesting that a 
significant part of the CME mass indeed comes from 
the dimming region. The mass flux carried by the out- 
flows has also been estimated from observations. 

Several temperature-dependent outflows have been 
found immediately outside the (deepest) dimming re- 
gions. The speed increases with temperature and it 
reaches ~150 km s _1 at log (77K)=6.3. Interestingly, 
our RB asymmetry analysis is able to detect some of 
these temperature-dependent outflows. These outflows 
are interpreted as evaporation flows which are perhaps 
driven by enhanced thermal conduction or nonthermal 
electron beams along reconnecting field lines, or in- 
duced by the interaction between the opened field lines 
in the dimming region and the closed loops in the sur- 
rounding plage region. 

Profiles of emission lines formed at coronal and 
transition region temperatures clearly exhibit two well- 
separated components in erupted CME loops and EUV 
jets. Besides an almost stationary component account- 
ing for the background emission, there is a highly 
blueshifted (~200 km s _1 ) component representing 
emission from the erupting material. The two com- 
ponents can be easily decomposed through a dou- 
ble Gaussian fit and we have diagnosed the electron 
density, performed a DEM analysis, and estimated 
the mass of the ejecta. Different properties of the 
two components suggest the importance of separat- 
ing emission from different sources when studying dy- 
namic events. Combining the speed of the blueshifted 



component and the projected speed of the ejecta from 
simultaneous imaging observations, we have calcu- 
lated the real speeds of the erupted CME loop and 
EUV jet. 
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and J.-S. He are supported by the National Natural 
Science Foundation of China (NSFC) under contracts 
40974105 and 41174148, respectively. H. Tian is sup- 
ported by the ASP Postdoctoral Fellowship Program 
of NCAR. The National Center for Atmospheric Re- 
search is sponsored by the National Science Founda- 
tion. 
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Table 2: Mass losses estimated from different methods and outflow densities for three dimming regions 



Obs. 


Scanning Period 


I st Method 


2 nd Method 


CME 


Outflow density / 


ID 




(g) 


(g) 


mass (g) 


log (N e /cm- 3 ) 


1 


2006 Dec 14 19:20-21:35 


1.9xl0 15 


8.0xl0 14 


3.6xl0 15 


7.0 


3 


2006 Dec 13 01:12-05:41 


4.1xl0 i5 


1.4xl0 15 


7.0xl0 i5 


7.1 


4 


2011 Jun 21 02:51-02:56 


l.lxlO 15 


5.0xl0 14 




6.8 
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A. Profile asymmetries not caused by blends or noise 

The high-speed outflowing component discussed in Section 3.1 is usually much weaker than the primary compo- 
nent, which may prompt people to think whether these high-speed outflows are caused by a weak blend or simply 
random noise. However, from the red contours in Figure Q] we can clearly see that most Fe XIII 202. 04A profiles 
showing significant blueward asymmetries actually have a high enough signal to noise ratio. Moreover, the blueward 
asymmetries form patches and they coincide with patches of significant blue shift and enhanced line width. If the 
blueward asymmetries are cause by random noise, we should see a random distribution of asymmetric line profiles in 
space. 

The cause of the blueward asymmetries by possible blends at the blue wings of the line profiles can also be ruled 
out, since we see these blueward asymmetries in not only one line, but all strong coronal lines in the EIS spectrum. 
As discussed in Section 3.1, the significant blueward asymmetries in the velocity range of 70-130 km s _1 at wings of 
these line profiles can not be caused by identified blends. 

There is another way to verify the above argument: the center to limb variation. We could not find any EIS obser- 
vation of a coronal dimming region as it rotates from the limb to disk center. However, we do have EIS observations 
of ARs as they rotate from the limb to disk center. As discussed in Section 3.1, the high-speed outflows in dimming 
regions seem to be very similar to those in AR boundaries. We take three E IS raster scans of th e AR 10978 from 2007 
Dec 10 to 15 for an analysis. Details of these observations can be found in lBryans et al.1 ( 2010|) . The AR was close to 



the east limb, disk center and west limb on Dec 10, 12 and 15, respectively. Figure [20] shows the spatial distributions 
of the single Gaussian parameters and profile asymmetries (averaged over the velocity interval of 70-130 km s _1 , as 
obtained from the RB p profiles) for the three scans. Center to limb variations of the profile asymmetries are clearly 
present. When the AR was on the disk center, we see prominent blueward asymmetries at both boundaries. As the AR 
rotated to the west limb, the profile asymmetries disappeared at the western boundary. And the profile asymmetries 
almost disappeared at the eastern boundary when the AR was close to the east limb. Clearly, it is hard to explain this 
phenomenon by noise or blends. While our scenario, namely a high-speed outflow superimposed on a nearly static 
coronal background, can easily explain this center to limb variation by taking into account the LOS projection effect. 
When the AR was close to the west limb, the magnetic field lines in the western boundary are almost perpendicular 
to the LOS so that the projection of the outflow speed on the LOS direction is very small, leading to a very small 
velocity offset between different components and greatly reduced blueward asymmetries in the velocity interval of 
70-130 km s _1 . This effect would also reduce the magnitude of the blue shift and line width of the total emission, as 
also revealed by Figure l20l 
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